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HATICaiAL ADVISORY CCMMITTEB FCE ifflROMDTICS 


TECHNICAL NOTE 2^97 


GEMRALIZEH) COKICAL-^TCW FIELDS IN 
SUPERSONIC WING THEORY 


By Harvard Lomax and Max. A. Heaslet 


SUMMAE?Y 


Linearized, ccmpresslBle— flow analysis is applied to the study of 
q.uasi-conical supersonic wing theory. Sin^e— integral equations are 
derived which relate either the loading to the shape of a lifting 
surface or the thickness of a symmetrical wing to the pressure dlstrihu:— 
tion for triangular wings with subsonic leading edges. The forma of 
these eq.uations and their inversions are simplified through the intro- 
duction of the finite part and the generalized principal part of an 
integral. 

Applications of the theory, in the lifting case, include previously 
known results. In the nonlifting case, it is shown that for a specified 
pressure distribution the theory does not always predict a unique thick- 
ness distribution. This is demonstrated for a triangular plan form hav- 
ing a constant pressure gradient in the stream direction. 


IHTRCfDDCTION. 


If a sufficiently thin wing at a small angle of attack is placed in 
a uniform stream, its aerodynamic properties can be determined by means 
of the analysis associated with linearized compressible— flow theory. 

If, moreover, a Cartesian coordinate system is used such that the wing 
is situated on or in the immediate vicinity of the xy plane and the 
stream flows parallel to and in the direction of the positive z axis, 
it follows that the basic eq.uation for the perturbation potential 
cp(x,y,z) can bo written in the form 

~ ^yy ~ *^zz ~ ® 

where = Mq^-I, Mq being the free— stream Mach number. 

The application of equation (l) to wing theory is essentially a 
mathematical problem involving the solution of a differential equation 
with given boundary conditions. Consistent with the assunq)tions of 
linearized compressible flow theory, or small-perturbation theory, the 
boundary conditions expressing the prescribed physical conditions are 
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glY0H always at z = 0 and.j as a conssqusncej louadary cond-itloiis as 
well as the solutions are superposable . 

The techniques used in the solution of wing problems are, for the 
most part, adaptations of existing mathematical methods to the specific 
type of boundary values and their suppcartlng surfaces that occur in 
aerodynamics. In particular, it is often possible in theory to mak e a 
reduction in the number of independent variables by virtue of known 
geometric or physical conditions. The conical— flow— field ana lysis of 
Busemann (reference l) provides in this way a means of descending from 
a three— to a two-dimensional potential equation. 

A conical flow field is one in which the perturbation velocity comr- 
ponents and the induced pressures are constant in. magnitude along any 
ray from the ax>sx of the field. In this case , the perturbation poten- 
tial may be written in the form 

)= ( 2 ) 

where cp is a homogeneous function of degree one in the three variables. 
An obvious generalization of this concept leads to the consideration of 
homogpneous potential fields of higher degree or, as they are eometlmes 
called, quasi-conical fields. If cp is homogeneous of degree K + 1, 
it follows that 


'P(3:,y,z) 



Eqmtion (2) yields coniceJ. velcxjlty fields, the degree of homiogenelty 
being zero, while for equation (3) the quasi-conical velocity fields are 
homogeneous and of degree k. Applications of these quasi-conical 
fields to pitching and rolling triangular wings have been given by 
Brown and Adams (reference 2), while Eibner (reference 3) has us§d sim- 
ilar methods in the consideration of cancellation elements. Further 
examples may be found in the literature. 

If new variables are introduced in equation (l) such that 


■^=11, ^ ^(z^y^z) = cx*''*'^ fl(TJj^) (^) 

where c is an arbitrary constant, the transformed partial differential 
equation is 

(ti^-I) iinq + 2 ti^ + (^^-1) - k[2ti + 2^ - (k+1) «]= 0 (5) 
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Thus, for a quasl-conical flow field with apex at the origin of the 
coordinate system, the resultant differential equation is elliptic for 
all •values of satisfying the inequality 

that is, for all points inside the foremost Mach cone 

X2 _ p2y£ _ pSgS 0 

The analysis of particular problems is therefore intlma-tely associa'ted 
■with the study of two-dimensional, elliptic— type equations and is espe- 
cially suited to the use of complex— ■variable theory. This is the 
approach taken by many investigators. In references 4 and Lagers-fcrom 
and Gennain have developed these methods in considerable de'bail. 

A different approach to the 
study of lifting surfaces in conical 
flow fields has been given by Brown 
(reference 6) and in reference 7. 

In this approach a basic lifting ele- 
ment carrying a uniform load dis-tn*i— 
but ion and extending radially from 
the apex of the field (see sketch) is 
considered first. The induced veloc- 
ity field is calculated for such an 
element lying in the plane of the 
wing and inclined to the stream 
direction at an arbi'brary angle 6, 

The solution of a particular problem 
then proceeds along one of t-wo lines. 

If the loading is given, the s-trength 
of each element is fixed and the cal- 
culation of the lifting-surface ^cme try .depends only upon carrying out 
the in'tegration. This is referred to as a direct problem. On the other 
hand, if the geome'fcry of the -wing is given, the loading is xtnknown and 
the strength of the elenKnts must be adjus-ted so that the resultant vei*— 
ticeil Induced velocities are consistent with the given wing slope every- 
where on the plan form. The solution of such a problem depends upon the 
Inversion of a relatively simple singular integral equation nnfl is 
referred to as an inverse problem. 

As ■was poin^fced out in reference 8, similar methods apply to non- 
lifting problems in conical flow fields and pressure did^tributions cor- 
responding to conical elements of thickness can be calcula^ted. In such 
cases, however, the direct problem, that is, the one involving the e^val— 
nation of an integral, is the one in which the slope of the wing surface 
is given; and the inverse problem, that is, the one involving the solu- 
tion of an integral equation is the one in which the shape of the pres- 
sure distribution is prescribed. 




/ 


/ 


/ 
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The present paper is concemed with the generalization of the 
basic elements of references 7 an^ 8 and their application to lifting 
and nonlifting problems in quasi— conical fields' of flow for cases 
involTing subsonic— type leading edg^a. Only solutions to the inverse 
problems -will "be considered and at all times these ■will be ob-bained by 
inverting the in'tsegral equation. 

The orders of singularities that arise in the analysis are such 
that it is convenient to use the concepts of the fini'te part and gener- 
alized principal part of improper in-tegrals. These ^neralizations will 
prove to be of importance for their ncbational efficiency and permit a 
simplified trea'fcment of the derivatives of singular integrals. Hadamard 
(reference 9) has pointed out clearly the necessary steps in the trea'b— 
HBnt of such improper integrals, but did not stress the role of the dif- 
ferential opemtion in ob'taining his in'tegrals. Since some differences 
exist be 'tween Hadamard* s definition of the finite part and the one used 
here, when extensions to multiple in'tegrals appear, a different no'bation, 
consistent with reference 10 has been adopted. The generalization of 
the principal part has also been discussed in reference 3» 


LIST OF IMPaRTAHT SYMBOLS 


ao 

0 ( 0 ) 

m 


Mo 

P 


Po 

P 

q. 


speed of sound in free stream 

load dis-bribution on lifting surface as a function of 0 
slope of radial element relative to free— stream direction 
slope of rig^t leading edge relative to free-stream direction 
slope of left leading edge relati-ve to free-s-tream direction 


free— stream Mach number 



local s-batic pressure 

free— stream s-batic pressure 

angular rate of roll in radians per second 


free— s-bream dynamic pressxire 



angular ra'be of pi-bch about wing vertex in radians per 
second 


Q 
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11 

V 


2 :,y,z 

a 

3 

6 


5 l 

Al 

Aw 

Ap/q 


free— stream Telocity 
streamwise perturtiatlon Telocity 
perturliation Telocity norml to plane of wing 
Cartesian coordinates introduced in equation ( 1 ) 
wing an^e of attack 
_ 1 


angle "between free— stream direction and line througji wing 
Tertez 

angle "between ri^t leading edge and stream direction 
angle "between left leading edge and stream direction 
discontinuity in u in plane’ of wing — U2) 

discontinuity in w in plane of wing (w^ ~ "W"^) 


load coefficient 


Cpz - Pu) 
1 


K 


X 

n 

Po 

T 

9U^yjz) 


constant determining degree of homogeneity of quasi-conical 
Telocity field 
(See eq.uation (3).) 

conical Taria"bles Introduced in eq,uation ( 4 ) 

m^, iiioP, miP 

slope of wing surface relatlTe to free— stream direction 
Mach angle (arc cot p) 
free-etream density 

region of integration in eq.uations (lO) (27) 

pertur"bation Telocity potential introduced in eq.uation (l) 

function related in eq.uation (4) to pertur"bation Telocity 
potential of a quasi-conical flow field 




6 


NACA TN 2497 


Subscripts 


u denotes conditions on upper stirface of wing 

Z denotes conditions on lower surface of wing 

AHALYSIS 


This investigation is confined to a consideration of inverse prob- 
lems; tbat isj problems that require tbe inversion of an Integreil eq.ua— 
tlon. As has been pointed out, these problems correspond to the two 

following cases: either the load 

distribution over a given lifting 
surface is to be determined or the 
thickness distribution corresponding 
to a prescribed pressure distribu- 
tion is to be calculated. The given 
conditions must, of course, be such 
that a q.ua8l— conical flow results. 
First, therefore, the plan form will 
be chosen, as shown in the accom- 
panying sketch, so as to have an 
apex at the origin of coordinates 
and to be of semi— infinite extent. 

The traces of the foremost Mach cone 
are inclined to the positive x axis 
at the Mach angle ±n = ±arc cot 3 
and, since only subsonic leading edges are being considered, the leading 
edges of the plan form are inclined at angles smaller in magnitude 
than 11 . Denoting these angles by 6© and Bi and measuring them from 
the X axis positively in the conventional counterclockwise direction, 
it follows that the eq.uations of the leading edges are 

y = X tan 8 q = m^x and y = x tan Bx = mxx (6) 

In the sketch Bq is positive while Bx is negative. 

The boundary conditions for the two types of problems may be stated 
as follows, where subscripts u and Z are used to denote conditions 
at z = Of z = 0-, respectively; 

Lifting case : Over all the xy plane Aw = w^^ — w^ = 0 and, 

except for the region occupied by the plan form, Au = Uy^ — U 2 = 0. 

On the plan form, vertical induced velocity is specified in either 
of the forms 
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''u = = 5^0 7^ ’’1 (f) (j) 

where 7i or 7z are known polynomials in z/y» 

g onlifting case : Over all the xy plane Zhi = = 0 

eind, except far the region occupied hy the plan form, 

/2iw = Wu — W2=0, On the plan form, the streamwise induced veloc- 
ity is specified in either of the forms 

= U2 = Vo Vi ^ = Vo x^ Vg 

where or are known polynomials in x/y. 

The solution to the two prohlems wil l he attained after con- 
sidering first a lifting element and a thickness element and then, 
for each of these prohlems, the basic integral equation is deter- 
mined hy smmning the appropriate elements. These derivations are 
given in the following sections. In small-perturhation theory the 
local load in coefficient form is related to u hy the expres- 
sions 


£p ^ PZ~Pu ^ ^ 
1 - V. 




(7) 


and the slope in the streamwise direction of an arhitrary surface 
z = z(x,y) is related to vertical Induced velocity hy the expres- 
sions 


Xu - 


^Zu _ 




hx V„ 




> 



(8) 


Lifting Case 


Upwash field of lifting element .— Consider a radial element emanat- 
ing from the origin and assume that the load carried hy the element is 
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( 9 ) 


Mo 


/ 

/ 

/ 

/ 

/ 

/ 

/ 



where C is a constant for a fixed 
position of the element. If the 
eleaent is inclined to the x axis 
at an angle 5, its upwash field can 
he calculated hy subtracting the 
induced fields of two triangular 
plan forms with vertex singles equal 
to 6 + A6 and each triangle 
having one side fixed for conven- 
ience along the x axis (see sketch) 
Assume, first, that B is positive. 
As shown in reference 10, the upwash 
field of the triangle can he ceilcu— 
lated frcm the known load dlstrlhu— 
tion hy maans of the fundamental 
formula 



(x-xi)^ (xi,yi) dxi 


kl~lxf J (x-xi)^-P^(y-yi)^ 


( 10 ) 


where the region t is the area on the plan form that lies ahead of 
the traces of the Ifech forecone from the point (x,y,0). The bars on the 
Integral sign indicate that the generalized principal part of the inte- 
gral is to he evaluated (see appendix). By definition, if 

r _ G-(y,yi) + constant (lla) 

o' (yi-y) 


is a known indefinite integral, the definite integral is evaluated as 
follows 


f = 0 ( 7 ^ 13 ) - G(y,a), a,h 4 7 (ill*) 

In case the singularity in the integrand lies outside the region of 
integration, the definition yields, of course, the , conventional definite 
integral. In- some cases in the following analysis the principal— part 
sign will he used to express relations valid for singularities both 
inside and outside the range of integration. 
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The Integration with respect to in equation (lO) is carried 

out over the area hounded hy the lines 

yi =0, xi = yi cot 5 = yi/m, (z-xi) = ±p(y-yx) 

andj with the suhstitutions, 

e = pm, Tj = py/z. Til = pyi/z (12) 

leads to the results 


1 wu pc r 


e(^+n) 

1+0 


(n-ni) ® 


O<0, -l<Tj<0 
(13a) 


4ir 


£ 


^(i-q) 

1-0 


'Hl*'^ /is Illv2 i . 

_ (r^^i) dTii; O<0, 


0 <TJ<1 
(I3t) 


If 0 and m are negative, the limits in equations (l3a) and (l3h) are 
reversed so that as the limits are now written changes in sign are 
required in the equations. 

In order to obtain the upwash field of the required lifting element, 
it is sufficient to perform a direct differentiation for, if of a 
plan form with vertex angle 6 is of the form 

w^ = f (0,Tj) 

it foUowB that when O<0, the value of upwash induced hy the element 
may he denoted dw^ and is 

dw^ = f(0+/50,ii) - f(0,q) = <10 

and when 0 < 0 


dwu = f(0,Ti) - f(0+^,q) = - g d0 


If this process is carried out in equations (I3a) add (l3h) and if the 
transformation 


41 
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is introduced, the eagpresslons for dw^ hecome 


\^J Vo 


/p f ^ _ pc de 

\xj Vo 4« 



- - - . -i<Ti<e<i 

(J:!)- - i< a< .< : 

(e-n ) (e-t t^ v^i-t2 


(l4a) 




Equations (l4a) and (l4h) provide the upwash fields for any radial ele- 
ment, regardless of the sign of 0. Integration by parts leads to the 
alternative forms 




dvg ^ _ pc d9(K+l) e*^ 

Vo 


dv^i _ _ pc da(K+l) 0 *^ 

Vo " 



(t]— t )*^ a/ 1— t^ dt 


(t)— t ) vi— t^ dt 
, ,K+2 ^ 

(0— t ) t 


— 1< •q< 0<1 


- 1<0<T)<1 


(15a) 


(I5t) 


vhere in one of the integrals the singularity in t req.uires the use of 
a Cauchy principal part. 

Derivation and inversion of Integral equation .— If now the lifting 
elements cover the region between 0x and 0o and C is a function 
of 0 determining the lift carried along the radial element at that 
point, upwash produced by the resultant plan form is given by the rela- 
tion 



P(k+1) 

4ir 


r ^ n' 

/ 0'' C(0) d0 / 

L Joi 


(a) aa r” at ] 

J-1 t J 


(n-t)'^A-t^ dt 
(0-t)'^"'^ t 


(l6a) 


This result can be written as 


/p^ ^ _ p(k+i) 

Vo 4« Jq^ 


e*^ C(0) H(0,ti) d0 


(l6b) 



MCA TN 2497 


11 


where 


H(0,71) 


7_a (e_tf+^ t 


dt. 


dt. 


9< 11 




The function H(0 jTi) has a sln^ile pole at 0 = i^, and the integral 
expression for w^ in equation (l6h) is therefore evaluated as a Cauchy 
principal part. 


The boundary condition to be satisfied by equations (l6a) and (l6b) 
is that ■Wu/x^ Yq Is a polynomial of degree k in the variable q. 
It follows that the (k+ 1) derivative of the ri^t-hand member of equa- 
tion (l6b) must vanish. 


Thus 


0 = ° <10 (IT) 

where use is made of the generalized principal part of an integral 
defined as (see appendix) 

A(xi) dxi ^ J:. A(xi) dxi 

Ja (xi-x)^+^ nl 


A(xa)lnlx-xi|dzi 

Here again the definition applies regardless of. the valiie of x but is 
of particular significance when x lies within the region of integra- 
tion. 


Continuing the calculation in equation (l7)j one has, after taking 
the derivatives with respect to q, 

0 = 0*^ G(9) id 

n ^01 
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wiilcli may "be written in the form 


f 9*^ c(e) de ^ Q 

\^/ v7ei 

The function C(0) is thus to he found through the inversion of the 
integral equation 


/' 

JQx 


O 0 *^ c(e) d9 

6-^ 


K 

i=o 


ai 


(18) 


The inversion of the Integral equation 


is known to he 


f(x) 



g(l) dg 
x-l ^ 


a< x< h 


g(x) 


1 

tj (h-z)(x-a) 



f(l)7(t-i)(^~a) 

x-i 



(19) 


where A is an arhitrary constant to he determined from physical con- 
siderations. ThuBj the solution to equation (l8) for 0i< 0< 0^ is 


0'* C(0) 


1 

-/(0Q-0)(0-0l) 



and this leads to the expression 


0*' C(0) = y (20) 

* </ (0o“9 ) (0-0 1) 

where the coefficients hi are functions of 9q and 0i hut not func- 
tions of 0. 

Belation of gene^ral solution to wing ge cone try- .— Brom equation (9) 
the loading on the plan form is, since 0 = Py/x, 
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Ap 


a 


K+l 


L=0 


'y(0Q~^ ) (0-0 x) 


( 21 ) 


wliere the coefficients hi must he determined ff’cm known Information 
about the surface geometry. 


Consider next the identity 





(h-t)(t-a) 



Ja (0— t) V(h-t)(t-a) 


0, 


a<0<h 


( 22 ) 


where i is zero or a poaitiTe integer. This expression implies the 
equality 



(q— t)*^A/l— t^ 
(0— t ) t 


dt = 0 


since the latter form can he broken, by expansion into rational frac- 
tions, into integrals like the left nsiaber of equation (22). The equal — 
ity 


J-i t 



(ri-t)*^yi-t^' 


dt 


follows where the principal-^art sign is needed on hut one side of the 
equation, depending on the valiie of 0 relative to q. K?om this 
result and equations (l6) the following relations are obtained: 




C(0) 



(q— t)*' Vl— t^ 

~ ,K+2 

t(0-t) 


dt. 


q<0 


(23a) 



Vo 




(q— t) V^l— t^ 
t(0-t)’^^ 


dt. 


q>0 


(23h) 
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T3ae range of tj in these two equa- 
tions has been restricted, respec- 
tively, to negative and positive 
values in order to avoid mathemat- 
ical difficulties arising when sin- 
gularities occur simultaneously at 
t = 0 and t = 0, The shaded por- 
tion of the accompanying sketch is 
the region of integration in equa- 
tion (23a). In the sketch, the 
inequalities — l^Si^O, O<0 q ^1 
have been assumed. This implies 
that the plan form has two subsonic- 
type leading edges and such a con- 
dition will be asBxaued to apply 
henceforth. 

It is particularly convenient 
to invert the order of integration 
in equations (23a) and (23b). 
in the forms 


f (n-t ^ ito 1 

J-l t(0-t)’''^ (K-fl)l \3a/ t(0-t-<x) 


(ri-t)*‘*/l-t^ _ iim 1 / ^ 1 --^ 

t(0-t)'^+2 ■ 0«+l)I VSa/ t(0-t-<t) 


substituting from equation (20), and Inverting the order of integration, 
leads one to the f onns 



11m -P ^ (T]-t)*^ Vl-t^ 



tt+i 

^ hi**©! d0 

i=o 

(Q— ■fc—otr) 


T1<0 



oalH 
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:5 


7o 


11m -P ( ^ (Ti-t)'^V^l-t^ 

a— >0 4jtKl \3a/ ^ t 


K+l 

) b. 0^ d0 

p^oh 

Joi (e-t-^)J (6^-9) (e-ex) ' 


[ 


11 >0 


Use of the algebra ic identity 


(9 


JH- = y (t+a)*’-'' + 

-t-a) A (0-t-a) 


J=i 


permits the expression of Vyj^ in the altematiTe forms 

^ . vK+l 


/pN _ 11m L A ) 

\^J To “ 0.-^0 4irKl \aa/ 


x: 


(n-t )*^ «/i— t^ 


dt 


y h. y 0 ^“J(t+a)«J“^d 0 

A U k 

Joi v'(0o-0)(0-0i) 


/: 


(Ti-t)'^Vi-t^dt r®o ^ 

<1 0 


K+l 

^ hjL (t-Hx)' 


d0 


01 (e-tMx) 7(00-0) (0-0 1) 


IC+1 

y (t+a)^ d0 


(11-4)*^ 7l-t^dt p^o i=o 

A- „ ■*= '^01 (0-t-^i) 770^l0T(0liiy 


'01-a 
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/ vK K+1 

( i ) = 11m _ P ( ^ ) 

Vs/ Yq a->o 4irK* V^/ 


K+l i 

I “ z 




J, t Jo. 


9i a/ (9o— 0 ) (0—0 iT 




(i^-t)*^ ^/l-t^dt r®o 1=0 

"t a. f Q—M—n 


K+1 

^ (t+a)- 


de 


01 (0— t— a) J (0Q— 0 ) (0— 0i) 


K+1 


y li (t+a)^ 

(T]-t)*^ 'N/l-t^dt /^o 


d0 


'00^ t 


&x (0— t— cc) ^/(0Q“0y(0H9l) 


The double integrals occurring first in the ri^t-hand members of these 
equatims are of degree k in a and their deriTatlves consequently 
Yanish, MoreoeYer, from the identities 



d0 

(0-z) J(b-0)(0-a) 




— TC 


a/ (x-a)(i-^) 


, b<s 


0 , a<x<b 

7t 


A/(b-i)(a-z) 


, x<a 


(24) 


it follows that the last tenns in the rl^t members vanish nnii that the 
expressions for w^^ become 



lim. ~P ^ —) 

a-^o IjKj V 


K+1 

V(0o--t^)(0i-t-a) 
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^ = lim f 

\^/ 7q a,->okKl\^J 

n^o-^ ll(t+cp)^ 

^ (t+a.-«Q)(t+a-^i) 

If, as in the appendix, the finite part of an integral is 

A(xi)dxi _ 2^ A(xi)dxj. _ 

J & (xi-b)^"^^^® 1.3...(2i-l) (xi-b)^^^ 

2^ / A(xi)dxi 

1.3...(2i-l) Ndh^ (xi-b)^^^ 

the expressions for may he -written in the fctrm 



— - (ti— t)^ «/l~t^ / A'^ 

^o 4Ki Ai 


K+x 

I *>1*^ 


i=o 


A(0Q-t)(0x-t) 


dt, Tj <0 


(25a) 



K+l 

X.. 

i=o 


(t-0o)('*^i) 


dt, 11 >0 (25h) 


Equations (25a) and (25h) are the fundamental equations for a 
lifting surface^and, since has heen assumed known as a polyno— 

miail in T| = 3y/x, it remains merely to determine the unknown coeffi- 
cients hj^ hy equating coefficients of q on both sides of the 
expressions. In the form given, it api>ears that for K>0 the number 
of equations obtainable exceeds the number of unde-termined coefficients, 
Ho general theorems of de-tenninancy have as yet been es-tabllshed as to 
the uniqueness of the solutions, but applications to be made la-ter will 
Indicate the techniques involved in calculating specific examples. In 
the laterally symmetrical case, -where -6 i=9q, solutions are easier to 
determine and the fundamen-tal eqxiations are 
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r^o 

4 k 1 t 


K+l 



11 <0 


K+l 


j_n^o 

4k! j; t 



(26a) 


(261) 


Nonlifting Case 


A radial element emanating from the origin is to he constructed 
such that it has a quasi— conical thicloiess distrihution' 

\i = c/ 

where G is a constant and is the streamwise slope of the element 

as defined in equations (8), The derivation of the induced pressure 
field associated with the element follows closely the analysis in the 
lifting case. Thus, a triangular plan form is first considered where 
one side is parallel to the stream direction and with a vertex angle B. 
From reference 10, pressure coefficient can he "written in the form 




(x-xx) dxi 

[ (x-xx)^-P^ (y-yi)^]®''^ 


(27) 


where the region t is the area on the plan form that lies ahead of the 
traces of the forecone from the point (x, y, 0) and the integration 
with respect to xx involves the finite part. 

In the notation of eq,uations (13), the analogues to equations 
(l5a) and (I5t) are 



2 C dg 9*^(k+1) t(Ti-t)'^ dt 

J_x (g_t)''+2^/IIt^ 


2 C d9 e'^(K-t-l) t(q-t)'^ dt 

(0-t)'"+^-/53^ 


, -1<T1<0<1 


, — K0<tl<l 


(28a) 


(28h) 
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If C Is a function of 9 and the thickness elenents cover the region 
between 9i and 0 q, pressure coefficient on the plan form is given by 
the esp)ression 



lim ~^(*^+l) ^ 



t(n-t)*^ dt ^ 



0'^C(0)d0 


^ t(T]-t)*^ dt 


( 29 ) 


The boundary conditions require to be a polynomial of 

degree k in tj. If the (K+l)st derivative of equation (29) is set 
equal to zero^ the relation 


0 = 




t(q— t)*' dt 
t('q— t)*^ dt 


holds, and after further differentiation reduces to 


0 = 



o 0*^C(0)d0 



1 r^o e^C{9)dJ9 
(K+1)]J0^ (0-^1 ) 


The function C(0) satisfies the same integral equation that arose in 
the lifting case. (See equation (l8).) The solution can therefore be 
■written, as in equation (21), in the form 


= C(0) y'^ 




ai0l 


( 30 ) 


The equivalence 


^ t(p-t)'^ dt 

J-X (0_t)'^+^'/l-t^ 



t(q-t)'^ dt 
(0-t)^+2-v/l-t2 


I 
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permits the revriting of equation ( 29 ) in the forms 



-2(k+1) 



t(T]— t)*' dt 


T] >0 


. ^ -2(k+1) 

'P Itp 




t(Tl-t)*^ dt 

(0_t)K+2Vl_t2 


T] < 0 


(31a) 


(31h) 


Suhstltution for 6*^0 (0) from equation (30), inversion of the order of 
integration, and use of eq.uations (2lf) leads to the fundauBntal rela- 
tions 




0 - ^ 

\^xj p 713?“ 


K+1 

I 

1=0 




*/ {9o ~^ ) (® ) 


K+l 


dt, T| <0 ( 32 a) 


G = 2 / o t ( i|--^ ^ / dt, q >0 (32h) 

V -y(t-0o)(t-0i) 

When 01 = -0o^ these equations hecome 




r 


\xy ^ pKi V dt / 


K+X 

i=o 


dt 


(33a) 


n - _2_ y ° ~fc(t-i)^ ( A 1 
\nj P pKjJ \ ^ 


K+l 


V t^-0o^ 


dt 


(331) 


The determination of the thickness distribution corresponding to a 
given pressure distribution can thus be obtained from the above equa- 
tions by equating coefficients of q and solving for the unknown 
coefficients aj^. Specific examples will serve to make the steps 
clearer; such i>roblems will be considered in the following section. 
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APPLICATIOHS 


Since homogeneous fields of low degree hare already received con— 
siderahle attentionj several results have heen puhlished previously. In 
the case of thickness prohlems corresponding to specified pressure dis— 
trihutionSj however j solutions have never, so far as is known, heen 
sought in terms of the given pressures. Bather, the thickness has heen 
assumed known and the resulting pressure distrihution calculated. This 
latter attack involves no question as to uniqueness and a one-to-one 
correspondence certainly exists. l*Jhen pressure is prescribed first, 
however, it hecomes necessary to consider the possihility of nonunique— 
ness. In two-dimensional, low— speed flow a freedom of choice is known 
to exist and leads to the introduction of purely circulatory flow which, 
in turn, provides the nechanism of lift. No analogue to this occurs in 
the lov^-speed, two-din^nsional, nonlifting case when the hody is smooth 
and is assumed to close. In the following developments a multiplicity 
of solutions will, however, occur in the nonlifting case and bodies with 
given pressure distributions retain a degree of freedom. 

Equations (25) and (26), together with equatlcn (21), suffice for 
the solution of quasi-conic^ lifting problems while equations (32) 
and (33), together with equation (30), apply to symmetrical wings. In 
the applications to follow, the division into lifting and nonlifting 
cases has been maintained. The detailed treatment of equations (26) and 
(33) can be further simplified if the problems are separated into cases 
involving symmetry and antisymmetry about the x axis of the Imposed 
boundary conditions. Suppose, first, that the given values of upwash 
rvr\(\ pressure coefficient in these equations are odd functions of t]. 

It follows from physical considerations that the loading or surface 
elope, respectively, will be an odd function of q and that conse- 
quently the unkn own coefficients ai or bi must vanish for even 
values of the subscript i. If the transformation t = — t is made in 
either equation (26b) or (33b), the pairs of equations (26a) and (26b) 
or (33a) and (33t) yield consistent sets of simultaneous linear equations 
that can be obtained from single equalities. Hence, for wu or Cp 
expressed in odd powers of q. 



or 
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vhere tlie notation [ k / z ] in the summation denotes tlie largest integer 
contained in k/ 2 . 


In tlie same fashion, a simplification can he achieved for symmetri- 
cal houndary conditions. Again using the hracfcet notation to Indicate 
the largest integer, the resulting equations become, when or Cp 

are expressed in even powers of tj. 



p^o (q— t )'^ 
4k! t 



[ Ck+i)/2] 

I 

J=o 






dt 


or 




C(k-m)/2] 

K-fl X 
t ^ JC*fl 

vat/ -/t2_0^2 

Wings With Load Distrihutions 


(36) 


(37) 


The yawed triwnpiil p-r -u-ing — IgquatjonH (21) and (25) lead directly 
to the determination of angle-of-attack loading on a yawed triangular 
wing. This solution is well known and was calculated in reference 7 hy 
a method which was a particular case of the present theory. 

The houndary cohditlons are that = — ToCt on the plan form, 
hence k = 0 and, from equation (21), 


_ ho+hx9 

^ -/(9o-0)(0-0i) 


( 38 ) 
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Eq.uation8 (25) lead to the equalities 


\/l— ^ bQ+hit 

t St ^/(t— 0 q) (t— 5 i) 


p p^o S 

t St 

4 J -I +. St. 


ho+tit 


4 ^-1 t St V(8o-t)(®i-t) 

After integration hy parts, these relations hecome 


Hence, 


where 


p r® 
p 


o (to+tit) dt 

t^ y (i-t^ ) ( t-80 ) 1 ) 

(hp+hit) dt 

t^ V (l-t^ ) (0o-t ) (01-t ) 


a = £ hp Lp(9pj Si) + t>i Li(Sp, 0i) 

4 L J 

~ ^ tp Lp(-0i, — 0p) - 1)1 Li(— 01 ^ -0o) 


I'O^^OJ ®i) ~ ~ J 


t^ v(l^^)(t^p)(t^i) 


Li(0p, 0i) = 


-r 


t V^(l-t2)(t-0o)(t-0i) 


It is apparent that equations ( 39 ) can he solved for hp and hi 
in terms of the functions Lp(0p, 0i) and Li(0p, 0i). Substitution of 
the values into equation ( 38 ) leads to the expression for load coeffi- 
cient 


^ ^ 2a / 2G 
q PE V 


_ (0p+0i) 0— 20p0i 

^ {9^—9 ) (0-0 1 ) - 
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wliere E Is th e complete elliptic integral of the secand kind vith 
modulus '/1-Q2 and 

g 1-9 q9i - ^/(l-eog)(l-9lg) [./(l-9o)(l+9i) - «/(l-9i)(l+eo) ] 

Oq~^i 2(0q-0i) 

In the particular case when 0i = -6 q, the value of G hecomes 0 q 
and the resultant loading on the unyawed wing hecomes 

Zip W 0Q® 

*1 PE a/^o2_02 

where the modulus of E is a/ 1-0^. This latter problem could, of 
course, have been solved directly from eq.uation (36) in a much simpler 
fashion. 

The rolling triangular wing .— Consider next the case of an unyawed 
triangular wing rolling about its axis of symmetry. If the angular rate 
of roll is P radians per second, the boundary conditions on the wing 
are that w^ = — Py. In this case k = 1 and, from eq.uation (21), the 
loading on the wing is given by the expression 


^ = S (42) 

where the coefficients bo and bs can be deleted since the loading 
must obviously be antisymmetrical. Since [k/s] is equal to zero in 
eq,uation (3^)j t>i is found to satisfy the relations 



(^3a) 


(43b) 


Performing one differentiation with respect to t, in each of the 
two integrals, and then integrating by parts, leads one to the expres- 
sions 


Vo 


Pbj.0Q^ r^o dt 

4 Ji t^-N/lI^ (t2-0Q2)®/^ 


Pbi0Q ^ r^o dt 

4 "4 4^ a/ (l-t^ ) (t^-0o^ ) 
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t dt 

The last equation reduces to an identity while the former one heconBS 
— = ^ ^ E \ _ phi 

Yq k 9q2 J 40^2 

where the ^ ulus of the complete elliptic integrals E and E is 
k .= Vl-9o^* aquations (42) and (44) the loading is 


2E-^o^ (E+e) 


(44) 


n^o 


0 - ti0o T 




^ = ...o ^ 


t 


1 


4p 9o^ X 9 



(45) 


and is in agreement with the results of reference 2. 

The pitching tr iangular wing .— If an unyawed triangular wing is 
pitching about its vertex, the boundary conditions became w„ = -Qz 
w^re Q is the rate of pitch in radians per second. Erom equa- 
tions (21) and (36), the loading and the relation involving the undeter- 
mined constants are 


^ _ X bp+bgg^ 

> q. P •J'q^-9'^ 

and 


ifQ ^ P (q-t) Vi-t^ C bg+bgt^ 

^O ^ ^ \ ^t / /v/t^— 0^^ 


since k = 1 and the loading must be symnBtrical about the x 
Equating coefficients of t] in the latter equation, one has 


axis. 





Vllb2 

\^y 


^o*^t^ dt 

(46a) 

l>o+liat2 

(46b) 



26 


MCA TN 2^97 


or, after carrying out one differentiation and integrating "by parte. 


0 = ta 






■/ 

1^1 


t dt 


0^ = -ba 


^ dt 

^r. 'Ja 


— (to+2ta0 


900 '-'1 


dt 


tV(l-t^)(t2-0o2) 


0^0 v'(l-t2)(t^_0/) 

(■bo+2ta0o^) ^ ^ 

900 4, 


4(l-t2)(t^-0o2) 


The first term in the first equation is zero and, after integration of 
the terms in the second equation^ the foUcwlng relatione are obtained 


0 = ho+ 20^2 ta 

= ho(M) + (E-0o^K)h2 

^o 

where the modulus of the complete elliptic integrals is k = ^/l-0^. 
The solutions of these simultaneous equations are 


^o = 


+8 0Q^h?Q 

Vo r 0o^K+ ( 1 - 200 ^ )E 1 


, ha = 


-h 






£ind the resultant loading on the wing is 


^ ^ 4x Q 2 0Q^ - 02 


m 


Differentially deflected triangular wing .— If the two sides of an 
unyaued triangular wing are deflected differentially, Tertical induced 
Telocity on the wing is 
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and loading nnist therefore he asymmetrical* It follows that K = 0 
and, from equations (21) and (3^)j 


^ _ hi 9 

a 


and 

®o hi t 

t ^ 

An Integration hy parte, in the latter expression, reduces the integral 
to a standard form. A further integration leads to the equality 



and load distrihution is therefore 


/ip _ 8a 9 q 0 


(49) 


Triangular wing with paraholic twist .— Consider, finally, an 
unyawed triangular wing twisted symmeirical.ly such that its Tertical 
induced Telocity is of the form 


^ = r ye 
Vo 

where r is a fixed constant. Since k = 2, the relations 

2 

^ ( xN hp+hg0e 


and 


ri]2 = 



(ij-t)^ / _S'\ 

t \ 


ho+het^ 

-^/te-0oe 


apply. Detailed analysis will he omitted in this case since it follows 
the same pattern of deTelopnent used in the earlier cases. Three, simul- 
taneous equations inyolTing the two uninowns hp and hg are obtained. 
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"but tlia equation relating coefficients of tj can be shown to Yanish 
identically and the two remaining equations yield a uniqiie answer. The 
ezpressiai for load distribution is 


a 


-C-t 

P VP/ 


[2eo^K-(i+9o^ )s]f [-(3-^0^ )K+H2~ef )E]e^ 


( 50 ) 




w here th e modulus of the ccanplete elliptic integrals K and E is 

Vl^o^. 


Wings With Thickness Distributions 


Triangular wing with uniform pressure .— In reference 11 , Squire 
couBldered certa.j-n tbic kne sH Ifltribnt.lnnH for synimetrlcal nonlifting 
wings in conical flow fields and calculated the resultant pressure dis- 
tribution. The first of Squire’s examples was a triangular plan form 
with a uniform pressure distribution. It is instructive to consider the 
inverse of this problem and to seek the wing assuming the pressure 
variation known. The plan form is symmetrically disposed with respect 
to the z axis or stream direction while the boundary conditions 
require that Cp Is a constant over the entire wing. In this case, 

K = 0 and from equations (30) and (37) the following relations hold 


^^u _ A _ 


ao 




t ao 


(51) 

(52) 


The value of Sq may be found in a manner quite similar to the one 
used in the previous examples. Thus, after the differentiation is per- 
formed, equation (52) becomes 

pCp _ P®o t^ dt ap ^ n^o t^ dt 


/(1-t^ ) (t^-0o^ ) 



MCA TN 2^97 


29 


After integration, this yields 


PCp _ ^ ^ ^ ao(g-S) 

2 9q b8o 1 - 00 ® 


(53) 


where the modiilTis of the elliptic integrals is vl— 0 q^. If Sq is 

eliminated from equations ( 51 ) and (53 )> the slope of the wing on the 
upper surface is 


cLzu _ pCp(l-^o^) 1 

2{K-E) ^/9o®-0® 


_ ^ (l— 0 Q® ) ^ 


2(Z-E) 


The ordinate of the upper surface results from the integration 



m 


and is 


2b? (K-^) 2^{K-E) 


m 


Triangular wing with linear pressure gradient #— It is now proposed 
to determine the thickness dlstrihution for an im^wed triangular wing 
for which pressure varies linearly in the streamwise direction. Set- 
ting Cp equal to hx, it follows from equation ( 30 ) that the slope of 
the upper surface is expressible as 


_ X ap+ag9^ 


since the solutim is obviously symmetrical about the stream axis. 
Since k = 1, equation (37) becomes 



ap+agt^ 


dt 


(56) 


Equation ( 56 ) is an identity and, eifter the coefficients of tj are 
equated, leads to the two relations 
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0 = — i— ( i-1 dt 

'-'x Vi-=fc^ \^y 0Q^, 

P^~b ^ / ^'N ap+agt^ 

2 Ji 


(57a) 

(57b) 


It is of interest to ccanpare the series of equations Just 
developed vith. the ccarrespcaidJjig equations in the prohlem of the pitchr- 
ing triangular ving. Foocniallyj the algebraic steps aie the same and it 
is to be e 2 i)ected that^ Just as in the case of equations (46), two 
simultaneous equations will be obtained and that their scdutions will 
provide the constants Sq and aa. In the present case, however, equa- 
tion (57a) can be shown to "vanish identically and aa a result only one 
equation in "tiwo unlmowna remains. This means that an infinite number of 
possible saLutions exists. The foil. owing calculations -will supply the 
necessary de"bails to confirm these remarks. 


Coneider equation (57a) and introduce the transformation 

"b^ = T 


= To 


Then, by means of the relations 




the equation becomes 

r,T 


0 = ao f 




ae + 


dT 

f2r 

a > 


Vl-T 



' »/r-To 

dr / 

^2t 


T 




Vt-Tq 

p'^o 

2T dT 


a /^To 


as 


2t, 




/ 


2T /^dT-J- r’'°/^dT 
V 1-T broJi V 1— T 

/ 


T dr 


^ (1-t) ( t-Tq) dro V (1-t)(t-To) J 


(58) 

(59) 
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Each of the ahove integrals can he evaluated directly and the hracketed 
terms are in both cases zero. Since similar integrals occur in problems 
of this type, however, it is warthwhlle to give. the following general 
formulas (n an Integer) 


r° j°dT — 

^1 ^/(1-t)('p-To) 


= 0.; = - _(14To), n=l; 


and 


It 


X 


(2n-l)l 
4 (n-1 ) Ini 


(To“ + 1) + 


n-x 


^ 2(2^l)l(2n-2J-l)I 




(n^j)IJI(3-l)l(i>-J-l)I ° 


, n>r (60a) 


y ^ dT = «, n = 0 ; = |( To® + 2To-3 ), n^l; 


rt(Tp-l) ■ (2n-l)lTp^ ^ (gB+l)I ^ 
2®°(n+l) L (n-i)lnl 2 n! n» 


n-i 

I 


J=i 


(23+l)I(2n-23-l)l ^ 

(n-j)I JI Jl(n-J-l)! ° 


n>l 


(60b) 


. It remains to calculate the terms in equation (57b). If the dif- 
ferential relations in eq,uation (59) used, the desired eapression is 


-3®b 




as 


^D3-5^Di + D2 

“ dTp 


(61) 


o 
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vliere 


Dx 


J>2 


= r'^°—==i=M= = - 2 r dn^ udu = I Tq^ - ^ (I+Tq^)® 

^1 Vt^t)(t-To) 


/^To ^ 

y;^ 


^ r 

yT(l-T)(T-To) 


AT ^ _2 


oLn^ udu = -2E 


Da = r^° ^ = -2 dn® udu = I; Tq Da + 4 (I+Tq) Di 

Jx A(1-t)(t-To) vVo 

anfl tlie modulus of the elliptic Integrals E and K is a/i— Tq. 
Direct calculatim gives for the coefficient of a^ the expression 


B (3-^q) - 2K 
(1-Tof 


and for the coefficient of as 

2 Tq E — Tq g — Tq^ K 

(X-Tof 

The value of h from equation (6l) thus estahlishes for pressure 
coefficient the value 


S „2 


2z 


p"(l-0o^f 


I ao 2g-E(3-flo^) + 


as 


(00^+00 ) K - 28 


o=^]} 


(62) 


where the modulus of the elliptic integrals is </l— 8^^ • 

gram equation (55), the ordinate of the upper surface of the wing 


is 


r 


dZu 


1 aoS^+aaP^y^ 


m m 


= — + — yf_^ (ao+2p^^8E ) arc cosh 

2P%2 2moSr ^ 


(63) 
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In reference 11, Squire conaicLered the thickness distrlhution that is 
ohtalned hy neglecting the arc hyperhoLic function in equation ( 63 )» 
His results correspond to the case when aQ is -Q^oc^as and are 
specifically 


= 


= 


aoz 


apx 




p2(i_ao2)2 


E(3-0o^) K - (4-28 2) u] 


(64a) 

(64h) 


If the wing is cut normal to the 
stream direction to form a trailing 
edge, a tr i a n gular plan form and an 
elliptical cross-section result as 
shown in the accompanying sketch^ 

If the root chord of the wing is 
Cq and the maximum thickness at the 
trailing edge is t, the constant 
ao In equation (64a) is equal to 
P^m^t/cc^ and the analytical 
expressions for the upper surface of 
the wing and pressure coefficient 
are 




It is apparent that a multiplicity of thickness distrihutions with 
the same pressure distrihutions must exist. Consider first the case 
when Uq in equations ( 62 ) and ( 63 ) is zero.- The surface ordinate wnrl 
pressure coefficient are, respectively. 


^ = 


“o 


arc cosh 


mox 


( 66 a) 
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C = 2agz 9 q^ 

P ' p2(l-0o^)2 

Tlie lateral section of such a vlng Is shown for particular veilues of 
Ho and p in the next slcetch (the curre denoted hj n = »), Along the 
root cho[rd the thiclmess is zero while the Tn^Ti-rniTtn thickness position 
occurs at the "value of mox/y satisfying equation 




(1+00^) K - 2E 


(66h) 


mox/y 


2 V(mox/y)^-l 


= arc cosh 


ffioX 


or at approximately 


moz/y = 1»3128 or y/md^ = 0,7617 

The pressure coefficients of equations (65b) and (66h) are iden- 
tical provided the equality 


ag _ iDpt [(3 ~Qq^) ^ — ( 4 — 20q^^] 

0 ^ 20q^Cq2 [( 1 + 00 ^) e _ 2E] 

holds. From equations (65a) and (66a) it follows that the surfaces 



(67a) 


(67b) 


have precisely the same pressure distrihution as does also any surface 
given by the relation 

2Z 2Z3_ + 2n2jg 
t t(l+n) 



where n is an arbitrary multiplicative factor 
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Cross sections In tlie x = Cq plane of sicrfaces given "by 
eq,uations ( 67 a) and. ( 67 b) are shown in the accompanying slcetch for 





Section BB Section A A 


Mq =^/2~ ^d. mo = 1 / 2 . Also includjed. are sections calculated, from 
equation ( 67 c) for n = 1 and. n = 1 / 5 . 

Since hod.ies with the san© pressure d.istrihution can he found., the 
possibility of comhining results and. getting a body ind.ucing no chan® 
in the. free-stream pressure should.be investigated* Itom equation (^), 
it follows that for Cp = 0 , the arbitrary constants and as must 

satisfy the relation 


^ = - gg - E (3-ep^ ) 
ao (1+00^ ) E - 2E] 


and, from eqmtion ( 63 ), the ordinate of the resulting surface is 
expressible in the form 






arc cosh 

J 


where 


P3 = P2/P1 
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and 

' Pi = (1+00^) ^ - 2E 

Pa = (3-9q^) ^ ® 

The surface produces a real vlng only If remains positive and this 
condition leads to the inequality 

(m^jx/y) yCm^x/y)^ - 1 / arc cosh m^jx/y ^ Pg 


The range of m<iE/y is from 1 to «> it is easy to show that the 

left memher has the lower limit 1 at mox/y equal to 1, On the 
other hand, Pg can he written in the form. 


1*3 


^ 

P2-2(1-9o2)(k _E) 


and, since K — E>0 for 0 q< 1, it appears lUcely that P^ is 
greater than one* A more detailed check shows, in fact, that Pg lies 
hetween 3 aTid 9 fpcm ■vdiich it follows that the inequality can never 
he satisfied in the neighhorhood of the leading edge and no real wing 
with zero pressure coefficient is possible. The variation of Pi, Pa, 
and Ps is shown in the acoon 5 )anylng slratch in which the variables are 
plotted as functiona of 
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CQNCI.IIDIBG EEMABES 


It has "been shown that the assuir^tlon of quasi— conical flow in a 
supersonic field transforms the basic partial differential equation for 
the perturbation potential to an elliptic— type equation in two inde- 
pendent -variables throu£^out the region inside the Mach cone» It is 
therefore not sui*prising that sdutiens of wing problems, for both the 
lifting and the nonlifting case, lead directly to the consideration of 
an in-begral equation (equation (l8)) of the type encountei*ed in -two- 
dimensional subsonic theory. In the analysis and the applications of 
this report, it is shown that. for a large class of specified conditions 
the known inversion of the in-tegral equation produces solutions that 
require s-trai^tforward in-tegrations and the sol-vlng of simul-taneous 
linear equations. 

An unusual feature of the resulting theory is the fact that a 
multiplicity of solutions may appear in a gi-ven problem. In re-fcrospect, 
this degree of freedom is not surprising since it is -well known that a 
null solution exists in -two-dimensional subsonic theory and appears in 
lifting problems in the form of a purely circiilatary flow. In the s-tudy 
of subsonic symnetrical profiles, this arbi-tai^iness in the solution 
occurs witen the gecaretry Of -the -wing is to be de-termlned from the dis— 
-fcribution of pressure exer-bed by the fluid. Since, however, closure of 
the -wing is necessary, an additional condition is gi-ven which estab- 
lishes uniqueness in much the same manner that the Hitta condition 
imposes uniqueness in the lifting case. In the consideration of super- 
sonic q-uasi-conical flow, similar conditions to de-fcermine uniqueness do 
not necessarily arise. Mathematically, the condition of tmiqueness is 
de-fcermined from the degree of dependence bet-ween straul-taneous linear 
equatians. 

One further remark concerning the analogy be-tween lifting and non^ 
lifting problems appears to be pertinent. In -fcwo-dimensional subsonic 
theory the in-fcegral relation between the perturbation velocities along 
the wing surface may be -writ-ben, in the lifting case, as 

3t X-Xl 

where c is chord length. In the symme-brical case, the interrelation 
is expressed in the form 


u(x,o) 



w(xi,o)dxi 


X— Xi 
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Tims, aside from, tlie factor p, a ccm5)lete 
mathematical emalysis of the two prohlems, 
associated with a flat plate at zero angle 


duality exists in the formal 
Hence, the circoilatory motion 
of attack (see sketch) is 





Lifting case Thickness case 

analogous to a slope distribution associated with zero pressure coeffi- 
cient, HoweYer, in three-dimensional supersonic wing theory the basic 
relations between w and u on the wing do not haTe the property of 
duality (see equations (lO) and (27)). The assun^tion of quasi— conical 
flow, howeYer, together with the restriction to triangular— type plan 
forms with subsonic leading edges, brings the study of lift and thick- 
ness into more general parallelism and a close similarity exists 
between the final expressions in equations (25) end (32). 


Ames Aeronautical Laboratory 

National AdTisory Committee for Aeronautics 
Moffett Field, Calif., July 10, 1951 
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APEEKDIX 

THE GSEKERALIZED PRINCIPAL PART AND FIHITE PART OB’ AN IHTEGEAL 


If the singular integrand of a convergent .improper Integral is dif- 
ferentiated formally, "without due regard for the singularity, the result- 
ing expression is, in ganeral, lD5)roper. In applied theory, however, 
the differentiation is usually to he performed upon the in"tegral itself 
and in this case a cai*eful trea"fcment of the entire expression leads to 
a finite ans"wer. The "bwo most common examples of such prohlems arise 
in the evaluation of the Cauchy principal part and Badamard’s f ini te 
part. The following de"7elopment indica"tes the manner in which these 
cases are extended to include multiple differentiations, (The gener- 
alized principal part shall he concerned "with in"be grands having singu- 
larities within the region of integration and of order n where n is 
a positive in"beger; the finite part, on the other hand, involves ln"be— 
grands with singularities at an end point of the region of in"tegration 
and of order n + l/2.) 

Consider first the evaluation of Cauchy’s principal part. In this 
case a single differentiation is used nnH the expression 

sM = ^ y A(xx) In |x-xil dxx, a<x<h (A1) 

becomes, for constants a and h, 

g(x) = f'’ AhilSi 

Ja 

Hiere the symhol on the in"tegral sign indica"tes that g(x) is to he 
evalua"fced hy a limiting process defined as follows 


g(x) 


lim 

e— 


r"- A(.,) ^ 


LJa 

'7x+e 2 : xx 


(A3) 


To assure the conver^nce of this integral it is sufficient hut not 
necessary to assume that A(x) is different iahle at the point xx = x 
and that elsewhere "within the legion of in"begration A(xx) is either 
continuous or possesses ln"fcegrahle singularities. The concept of the 
Cauchy principal part is so well known that the symbol on the integral 
is often omit"ted, as shall he done here. 
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Turning next to the case of multiple differentiations, consider the 
expressions 


I = 






A 

3x X 



a< x< h 


(A4) 


where the limitB of integration are Independent of x and the symbols 
on the two latter integrals indicate that the generalized principal paid 
is to be calculated. From the definition of I in the first integral, 
it follows that 


T _ S -j b— X _ b— a 

^ x-a (b— x) (z-a) 


(A5) 


The simple definition given by eq.uation (A4) can be generalized to 
include integrals of tiie type 


la = A f A(zi) In |z-xil dxi = A P 

3xf az xi-x 

> 

T A(xi) dxi 
(zi-x)^ 


(A6) 


Eq.-uation (A6) defines the symbol appearing in the final member. It is 
possible, however, to relate this integration to the particular Integrand 
in the final integral of eq.uation (A4) by writing I 2 in the form 



/ 


A(zi) -A(x) 

Xi— X 


dxi + A(x) 


£ 


dxix 

Zl— z 


Then if A(xx) is integrable and if, at Xx = x, its derivative exists 
and is single valued, the expression for Is becomes 


la 


A(xx) - A(x) 

A ( 3 : 1 - 3 :)^ 


dxx 


+ A(x) 


b 

f 

Ja {xi-xf 



A(zx) - A(x) 

(3:i-z)^ 


A(x) (b-a) 
(b-x) (x-a) 


(A?) 


where the results of eq.uation (A5) have been used. 
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Tile first Integral in eq_uation (A7) is now in a form, that involves 
no extension beyond the concept of Cauc^'s principal part the 
evaluation of I2 may be carried out with that form. Furthermore, it 
can be shown from equation (A7) that if the indefinite integral of 
A(xi)/(xi— x)^ exists such that 



A(xi)dxi 

(xi-x)2 


G(xi,x) + C 


(A8) 


then the value of I2 can be found by following the conventional rules 
for substitution of limits so that 



A(xi)dxi 

(xi-x)^ 


= G(b,x) - G(a,x) 


(A9) 


The extension to hl^er ordered derivatives Is obvious. Thus, for 
a and b Independent of x, one 


I 


nfi 


-r>fT J A(xi) 2n lx-Xi| dxi = — 
ds '-'a ^2;^ 

n» A(xi)dxi 



A(xi)dxi 

(xi-x) 


(AlO) 


Equation (AlO) defines the final member appearing in it. The quantity 
can also be written in the form 


“•n+i 




b , ^b 

A(xi) - B(x,xi) ^ f B(x,xi)dxi 


(xi-x) 


n+i 


a (xi-x) 


n+i 


(All) 


where 


B(x,xi) = A(x) t (xi-x) + . . . + — 


(1^1) f^\ Q-i 

, ( 3 :i-x) 


and 


1 b 

r\ 




dxi 


1 1 
+ (-1) ^ 


^ ^ L (b— x)^ (x-a) -I 


,1< n 
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lEhe first n derivatlTes of A(3 :x) assumed to exist and "be single 

■rolued at = z wMle elsewtiere in the reinge of integration A(z) 
m£Ly possess integrahle singularities. The generalization of equation (A9) 
holds so that if 


/ + C (A12) 

then 

■b 

J (A13) 

a (xi-z) 

It is also possible to extend the definition of equation (AlO) to 
include a functional dependency on x in the numerator of the integrands. 
ThuSj replacing A(xi) hy A(x,Xi), equation (AlO) again defines uniquely 
a principal part integral proTided the first n deriYatives of A(x,Xi) 
vith respect to x and Xi exist at Xi = x. 


The original concept of the finite part was uaed hy Hadamard in 
connection with square root singularities. Consider the expressions 






h 

1 S 

2 Ja (^Xi 


(Al4) 


From the first integral in this relation it follows that 

h 


X 


dxi 


(‘b-xi)°'^^ -y/h-a 


(A15) 


The natural extension of this idea is to consider 

r = A(xi)dxi _ ^ A(xi)dxi ^ _ 1 A(xi)dxi 

ah Ja ^ ^ ^ 

where A(xx) is continuous at Xx = h and is integrahle elsewhere in 
the range of integration. 

The eYaluation of Ja can he related more closely to the integral 
of equation (Al4) after rewriting Jz in the form 
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Ja = ^ 


A(xi) - A(bj^ ^ r 


(A17) 


It follows tliat 


(t>— Xx) ("b— Xx) ®- (t— s:i) 


(ai8) 


An interesting integration tecliniq.ue can "be evolTed from equation 
(Al8). Setting 

A(xx)dxx _ lim A(xx) -A(b) , _ 2A(1 i) 

and setting the indefinite integral of 


r A(xx)dxx 

equal to F(h,Xx) + C, it follows that 

■T = - CF(b,a) + C] 

'^a (l>-xx) ^ 

where 


lim r 2A(b) 

L Vl>-xx 


F(h,xx) 


ThuSj with the proper choice of the constant of integrationj the defi- 
nite integral is found hy substituting convent ionally the lower limit. 
In practice, C is often zero. 


Defining Jq+x ^ form 


‘^n+i 


■ aJ i 


■f 

a a/B-Xx a 


A(xx)dxx 




A(xx) 

a/S-Xx 


(gjQ— 1) A(xx)dxx 


2^ '4 (b-xx)^^^'^ 


(A19) 
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it follows that 



A(zi)dxi 



A(xi) - B(h,xi) 
(■b-xi)^"^ 



B(‘bjXi)dxi 

(l>-xi)“'*i' 


(A20) 


where 


B(h,xi) = A(h) - A*(h)(h-xi) + ... + ( _ ^ _ (h) , (i>_3.i)n i 

(rt-1)! 

and 

dxi ^ (-l)V / dxi ^ 1 

J& (h-xi)^'^l- • 1.3...(2i-l) Ja. 0>-Xi)l 21-1 (h-a)^*^. 


It is furthermore posslhle to extend the definition of eq.uation 
(AI 9 ) to Include a functional dependence on h of the Integrand A{xi). 
^placing A(xi) hj A(h,xi), equation (AI 9 ) again defines imiquely a 
finite part integral provided that 


11m ^A(h,xi) 


-v/b-xi 


= 0 
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